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Abstract 



We find a self-consistent pp-gravitational shock wave solution to the semi- 
classical Einstein equations resulting from the 1/N approach to the effective 
action. We model the renormalized matter stress-energy-momentum tensor 
by N massless scalar fields in the Minkowski vacuum plus a classical particle. 
We show that quantum effects generate a milder singularity at the position 
of the particle than the classical solution, but the singularity does not disap- 
pear. At large distances from the particle, the quantum correction decreases 
slowly, as 1/p 2 (p being the distance to the particle in the shock wave plane). 
We argue that this large distance correction is a necessary consequence of 
quantum gravity 
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I. INTRODUCTION 



The backreaction of quantum fields on the spacetime geometry can be described by the 
so called Semiclassical Einstein Equation (SEE) [0J 

+ aH^ + 01^ + A 9lu , = 8vrG (T;|f ss + (T^)) (1.1) 

where one includes as a source the mean value of the energy momentum tensor of the 
quantum fields. The higher derivative terms 

= —2R.^ V + 2g fMV OR — -g^ u R 2 + 2RR ilv 

V = -2R«„ a + URp, + ^g^UR + 2R«R av - ^g^R a pR aP (1.2) 

are needed to renormalize the theory. The dressed coupling constants G, A, a, (3 depend on 
an energy scale fi introduced by the regularization method. The dependence is given by the 
renormalization group equations. The mean value (T^) also depend on fi in such a way that 
the full equation is //-independent. The numerical values of a and (3 should be determined 
by experiments or could be computed from a more fundamental theory. They should be of 
the order of M 2 Planck for fi ~ M P[anck . 

Although a complete theory of quantum gravity from which the SEE could be derived is 
still lacking, there are several arguments that lead formally to this equation. On one hand, 
one can expand the metric around a classical solution and quantize both the linearized 
metric fluctuations and the matter fields on the classical background |TJ. The one loop 
effective action for this field theory gives the leading correction (in powers of h) to the 
classical Einstein equation. In this approximation, (T^) contains both graviton and matter 
contributions. As two loop corrections will induce additional terms proportional to ft 2 , in 
this context it is reasonable to look only for solutions to the SEE which are perturbative in 
ft, and therefore close to the classical solution. On the other hand, the SEE can be derived 
using a 1/N approximation. Assume there are A" decoupled quantum matter fields. Taking 
the limit N —>■ oo with GN fixed (and rescaling also the classical energy momentum tensor 
so that GTf^ ss is ^-independent for large N), one arrives @ at an equation of the form 
flLip. In this case, there is no graviton contribution to (T^), since it is suppressed by a 
factor 1/N with respect to the contribution of the A" matter fields. Using this point of 
view, the SEE is valid to all orders in ft, and it makes sense to look for exact, self-consistent 
solutions. However, it should be kept in mind that some of the self-consistent solutions to 
the SEE may be non-physical. It is not clear that one can ignore graviton effects at Planck 
curvature, even for large N. Moreover, Eq. (|1.1|) contains time derivatives of order higher 
than two, and therefore 'runaway' solutions could arise. We refer the reader to Ref. for 
discussions on this important issue. 

Physical, exact solutions to the SEE are interesting because they can show important 
deviations from the classical behavior. The standard lore about smearing of classical sin- 
gularities due to quantum effects should be proved in this context. The problem of the 
endpoint of Hawking radiation should, and in fact it is, being analyzed using the SEE, al- 
though progress has been important only in two dimensions ||. In four dimensions, (T^) 
is a complicated, non-local function of the geometry, and the SEE has been solved in a few 
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cosmological cases |J. In this paper we are going to find a plane fronted with parallel rays 
(pp) ~~ gravitational shock wave self-consistent solution to the SEE. From the exact solution 
we will be able to discuss the problem of the smearing of the classical singularity. We will 
also compare the exact and perturbative solutions. 

The paper is organized as follows. In Section II we introduce the pp-gravitational shock 
waves, and prove some useful geometric identities. In Section III we use a resummation of 
the Schwinger DeWitt expansion to write down the SEE for this particular metric. The SEE 
is solved in Section IV. We discuss our results in Section V and leave for the Appendix the 
case of the ultrarelativistic straight cosmic string. 



II. THE GRAVITATIONAL SHOCK WAVE 

It was stressed by 't Hooft || that at Planckian energies the relevant interaction among 
particles is the gravitational one, which takes a particular simple form. In fact the gravi- 
tational field generated by a highly energetic particle can be described by the metric of a 
gravitational shock wave 

ds 2 = —du dv + f(x±)S(u) du 2 + dx 2 ± , (2.1) 

where u = t — z and v = t + z are the usual null coordinates, and dx\ is the two-dimensional 
flat metric. The wave propagates on a flat background along the trajectory of the particle, 
i.e. u — 0. The "profile" function f(x±) completely characterizes the spacetime. 

Geodesies that cross the shock wave suffer a sudden shift in its v coordinate given by 

Av = —f(pi) , (2.2) 

where pi is the distance to the inplunging geodesic from the source on the surface u — 0. 
This effect is represented in Fig. [I]. 

Geodesies also have a refraction effect given by M 



cot(0 in ) + cot(9 re ) = -d p f( Pi ) , (2.3) 



where 6i n and 9 re are the incident and refracted angles respectively. This effect is depicted 
in Fig. @. 

The SEE ( |1 . 1| ) in principle involve terms containing powers of the curvature tensor, 
which in turn would be proportional to powers of the delta-function 5(u). These powers are 
meaningless in the context of classical distribution theory. In order to avoid this problem, 



we will consider a more general metric [13 



ds 2 = -du dv + H(x±, u) du 2 + dx\ . (2.4) 

We will show that for an arbitrary smooth function H(x±,u) the SEE becomes linear in 
the curvature and, at the end of the calculation, we will consider, the limit H n (x±,u) — > 
f(x±)5(u) ({H n } is a convergent sequence of well-behaved test functions). 

The curvature scalars of the metric ( |2.4j) formed out of the Ricci squared, Kretschmann, 
and Curvature scalar tensors, all identically vanish 
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R = , R a pR a ^ = , R a [3jsR a ^^ = 



(2.5) 



This does not mean, however, that the spacetime is free from any singularity. In fact, 
the non-vanishing components of the Riemann tensor are (apart from others obtained by 
symmetry properties) 

Riuju = --8^H(x±,u) , (2.6) 

which may still diverge (as we will next see). 

From Eq. (|2.6|) we see that the only non-vanishing component of the Ricci tensor is 

Run = ~\vlH . (2.7) 

Note that the only non vanishing components of the inverse metric are g uv ,g vu ,g n and 
g vv . Also that the Christoffel symbols with a subindex v,T% v , vanish. This implies that the 
covariant derivative V„ coincides with the ordinary derivative when applied to a covariant 
index of a tensor, and vanishes when acting on a u-independent function. 

From the above properties it is easy to show that any scalar formed out of the Riemann 
tensor and its derivatives must vanish. Indeed, a generic scalar can be built up as 

;/ ; ' : /^..//lV A T A ,..V A JUiiV.--V/?...)...(V...V/L.) . (2.8) 

The indices Tj can be contracted with any of the covariant indices of the tensors in 
(V... VR...). As the only non- vanishing components of the Riemann tensor are those given 
by Eq. ( |2.6j ), either \i or v should equal u. Therefore, T\ or r 2 should equal v. As R v ... = 
and V„-R... = 0, the scalar (|2.8|) must vanish. Note that what we have shown is that any 
expression containing Rp, vpa vanishes if both indices [iv (or, by symmetry, pa) are contracted 
with any other index. 

With similar arguments one can show that the most general tensor with two covariant 
indices must be a linear combination of g^ v and Fi^J^R^ [with F an arbitrary function]. 
The proof is as follows. Consider a generic tensor (2, 0) 

^^...^-(V Al V A2 ...V An i?^ (T )(V...Vi?...)...(V...V J R...) (2.9) 

The difference with respect to the scalar case is that now we can have up to two non- 
contracted indices in the Riemann tensor. If Rp, vpa has three or four indices contracted, the 
previous arguments imply that the tensor must vanish. Therefore, Rp Vpa must have exactly 
two non-contracted indices and, as we are considering a (2, 0) tensor, it must be linear in 
Rjivpu- Fo r r i and r 2 there are two alternatives: 

a) if T\ is contracted with r 2 , the tensor becomes V Al V \ 2 ...V \ n R ua . The indices Aj must 
be contracted among themselves, and one finds something proportional to (V\)2R ua . 

b) if T\ and r 2 are contracted with two derivatives Ai and A 2 , the result is proportional 
to "SJ\R VU . The other derivatives must be contracted among themselves, and one finds again 
the same result. 

The conclusion of this section is that any geometrical tensor built up from the metric, with 
two indices, and in particular, as we will soon see (Tp U ) (when the state is the Minkowski's 
in- vacuum state), must be a linear combination of g^ v and a function of acting on Rp V . 
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III. THE SEMICLASSICAL EINSTEIN EQUATION 



Let us write Eq. ( |1 . 1| ) for the particular case of the shock wave metric ( |2.1| ). Using the 



fact that H. 



u 1 ' 



0, /, 



hi/ 



R 



V 2 



R^ u , we find 



+ JV 2 R 



111' 



8ttG 



■ [11/ 



(3.1) 



Here (as in Eq. ( |1.1| )), we have split the stress-energy-momentum tensor into two pieces: a 
classical source plus the renormalized expectation value of the T^ u due to iV quantum matter 
fields. Since the state in which we take the expectation value is Minkowski's vacuum, these 
fields do not contribute to the classical source and, since the fluctuations of the classical 
source are suppressed by a factor 1/N (as are the graviton contributions) neither of these 
two terms appear in ( [3.1|) [For more details, see, for example, Refs. ||[|.] 

In general, we expect the function F to depend on the mass of the quantum fields and 
on the scale \x introduced by the regularization method. For massive fields, this function 
can be expanded in powers of (V^/m 2 ) (Schwinger-DeWitt expansion). For massless fields 
this expansion is inadequate, and F becomes a non-analytic function of Vj_. All this can be 
checked from the explicit calculations done by Vilkovisky and collaborators [TO] . They com- 
puted the effective action up to quadratic terms in the curvature, and the energy momentum 
tensor up to linear terms in the curvature. But for the pp-shock wave we have shown this 
is an exact result! The energy momentum tensor is linear in the curvature. Therefore, the 
function F is given by fll0,12| 



F(V 



m 



167T 



1. 
6 



+ 



V 2 



384tt 2 



dt t 4 In 



m 2 - |(1 



t 2 )V 2 



As anticipated, if m 2 ^ 0, it is possible to expand F in powers of (V^/m 2 



F(Vi) 



m 

167T 1 



m 



(f-«)MT5-] + 



6 



384tt 2 







2 " 

m 


f 1 dt t 4 


In 
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lv 2 \ 



the Schwinger DeWitt expansion is recovered. 

On the other hand, if m 2 = 0, the form factor reads 



1 V 2 



77T 



+ o((- 



(3.2) 



(3.3) 



F(Vi) 



1920tt 2 



In 



/// 



+ const V 



(3.4) 



and it is a non-analytic function of the Laplacian. 

It is interesting to note that in the massless case the form of the function F can be 
derived using general arguments ||11|| . Indeed, from dimensional analysis we know that F 
must be of the form F(fi 2 , Vj_) = V 2 L F(V 2 L //i 2 ). Moreover, G is independent of \i and (3 
depends on /x according to the renormalization group equation [|T^ 

dp 1 



V 

8ttG dfi " 960^^ 



(3.5) 
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As the whole semiclassical equation must be independent of /i, we must have pdF/dp = 
— l/(9607r 2 ). This equation fixes the function F up to a term that can be absorbed into a 
redefinition of f3, and the result coincides with Eq. (|3.4|). 

Now we are ready to write the differential equation for the profile function. Assuming a 
point-like classical source [] 



T, 



class 



pS(x ± )5(u) 



we obtain 



Vi+ [/3(fi)-Aln 




f{p) = 8nGp6(x ± ) ; A 



G 



240tt 



(3.6) 



(3.7) 



Note that the //-independence of the above equation is explicit: Eq . (|375| ) implies that the 
parenthesis multiplying in Eq. (|3.7|) is independent of p. Eq. (pTT|) has been solved in 
Ref. [|T3[ for the case ^4 = (i.e. for the fourth order theory). 



IV. SOLUTION TO THE SEMICLASSICAL EQUATION 



Plane gravitational shock waves (|2.1|) have the nice property of linearizing field equations. 
This feature allowed us to write Eq. (|3.7p in an explicit form. Again, this linearization allow 
us to seek for the solution to the field equation by using Fourier transform methods. The 
result is 



f(x±) = fan + 16vrGp 



dk 
4^2 



,,ik-x i 





'?- 






1 - 


(3k 2 


+ Ak 2 In 


m] 



(4.1) 



where Jgr = — 8Gp ln(p/p ) is the general relativistic result [14[], with p an arbitrary 
constant. The second term in the above equation contains the corrections to the classical 
result. 

To perform this integral we first note that since the problem has axial symmetry it is 
convenient to pass from Cartesian to polar coordinates in the perpendicular space, i.e. from 
xj_ to (p, if). Doing so, the angular integral in Eq. (|4.1| ) is a representation of the Bessel 
function fl5| J 



f(p) = f G R + SGp / dk kJ (kp) 





'(3-Aln(%) 




1 - 


f3k 2 + Ak 2 In 





(4.2) 



Note that in the fourth order theory (^4 = 0, (3^0), the integrand above has no real poles 
when the so called non-tachyon constraint, (3 < 0, is satisfied. When A ^ 0, the condition 
for absence of real poles is Ap, 2 exp((3 / A) < e. 



1 This is an approximation, but if one wants to deal with an extended source, our results below 
give the Green's function for that case. 
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In this integral (now on the magnitude of k) , if A 7^ we can absorb the term proportional 
to (3 into a redefinition of the scale p 



f(p) = fon-*Gp dk- 
Jo I 



A\n(k 2 fi~ 2 e-P/ A ) 
YAk 2 ln{k 2 fi- 2 e-P/ A ) 



kJ (kp) 



(4.3) 



Let us now discuss some properties of the self-consistent profile function. We first note 
that it is finite at the origin p = 0. Indeed, we can write 



f(p) = -8pG 

= fan + f Q + Af. 



ln(p/p ) + Mp/VA) - I dkJ (kp*) I - fc ^f^L 

11 \ 1 + k 2 In ( -^J 



k 



1 + k 2 



(4.4) 



where we have taken dimensionless variables k=\J~Ak, p*=pj \[A and A*=Ap? exp(fi/A) for 
A 7^ 0. Also, we identified with the subindices U GR\ the general relativistic result, with 
the corrections induced by terms quadratic in the curvature (proportional to K ), and 
with Af the higher order corrections. 

Using the asymptotic behavior of the modified Bessel function Kq and noting that the 
last integral in Eq. ( |4.4j ) is finite in the limit p — * one can easily check that f(p) is finite 
in this limit. This confirms the results of Ref. |13j where the first two addends of Eq. ( |4.4| ) 
have been obtained. 

The next question is about the singularity at the origin. Although all curvature scalars 
vanish for the pp-wave, the curvature tensor itself may diverge. In fact, 



TDV 

^pup 



S(u)d 2 f(p), R\ 



2p 



dpf(p) 



R p 

puu 



(4.5) 



that for the classical solution diverge as d 2 fcR ~ P^dpfcn ~ p~ 2 ■ When quadratic terms in 
the curvature are included in the classical action, the divergence is milder, the components 
of the curvature tensor diverge logarithmically |L3| d 2 (fcR + fo) ~ P _1 ^p(/gr + Iq) ~ l n P- 
We will now show that the divergence is still milder for the self-consistent solution. The 
second derivative of the profile function can be written as 



d 2 J = 8pG 



dk 



kJ'o'jkp*) 
1 + k 2 In J 



(4.6) 



where the primes denote derivatives with respect to the argument. Note that since J' Q ' (0) 
— 1/2, this integral diverges as the ln[ln()] of the upper limit. 
To see the dependence with p we split this integral as 



d''J ~ SpG 



1 + k 2 In J 



+ 



kln¥- 



h) 



A* 



(4.7) 



where k is such that k^ In 3> 1. The first term in Eq. ( |4.7| ) is finite in the limit p — >• 0. The 



second term gives, after introducing the new variable y 
by parts 



kp* and performing an integration 



7 



ik 



00 dy J%(y) 



In 



A*p 



1 9° 



^ J hop* 



In 



.4*p 



J"(fc p*) ln[ln(— - -)] - ln[ln(fc M*)] 



2" v " ur 7 l"" r "M*p* 2 ' 



(4.8) 



Using that J"(0) = —1/2, we see that <9 2 / ~ 2pGln[— ln(p*)] as p — > 0. A similar result can 
be shown for p~ l d p f . The conclusion is that the curvature diverges like ln[— In p*] near the 
position of the particle. 

What happens at large distances? In this case we expect the self-consistent solution to 
approach the solution perturbative in H (perturbative in A, in our notation). From Eq. ( |4.4|) 
we see that, to lowest order in A (for large p), 



Af(p) ~ -8 P G j o dk U (kp*) In ( — 



(4.9) 



This integral is ill defined and must be treated as a distribution. It is essentially the Fourier 
transform of ln(fc 2 /^4*) in two dimensions. To compute it, we introduce an integral repre- 
sentation for the logarithm 



Af(p) ~ 8G P / dz dk U (kp*) 
Jo Jo 



1 



1 



z + k 2 z + A* 



(4.10) 



The A*— dependent term gives a contribution proportional to S(x±) which is irrelevant at 
large distances. Omitting this term we obtain 



A/(p)~8Gp/ dzK {^P*) 
Jo 



!6GpA 



(4.11) 



Therefore, we expect that, at large distances, /(p) = faR + 16Gp^4/p 2 . It is important 
to stress that the perturbative and self-consistent solutions are completely different unless 
p — > 00. 

We have numerically computed the self-consistent solution starting from Eq. (|4.4j ), and 
confirmed the properties of the profile function both at the p — > and p — > 00 limits. The 
results for Af and d 2 p f are shown in Figs. (|3|) and (|4|) . In those figures we only show curves 
for A* < e, since Af diverges for A* = e. This is due to the fact that, as we already pointed 
out, the integrand in Eq. (|4.1|) develops a double real pole as ^4* approaches e, that becomes 
two single real poles for A* > e. These two real poles, k p , change the asymptotic behavior of 
Af for large p*, making it to decrease slowly (like a linear combination of Jo(k p p), p = 1,2). 
However, from the definition of A* right after Eq. ( |4.4j ), we and the renormalization group 
equation fl3.5|) , we see that A* is actually independent of p. Then, from the renormalized 
value of A, given in Eq. (|3.7f) , and the no tachyon constraint that implies (3 < 0, we can 
infer that A* < 1. 



V. DISCUSSION 

We have found a self-consistent pp-shock wave solution to the SEE in the 1/N approxi- 
mation. 
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We have shown that quantum effects make milder the singularity of the classical solution 
at the origin. In the classical theory, the curvature diverges like 1/p 2 . In the quadratic theory 
it diverges like ln(p 2 ). When the backreaction effects of massless quantum scalar fields are 
taken into account, it diverges like ln[— ln(p 2 )]. Therefore, the self-consistent solution can 
be trusted up to a distance very close to the position of the particle, where the curvature 
becomes of the order of Planck curvature. 

At large distances, the situation is similar to the case of the corrections to the Newto- 
nian potential: while classical terms quadratic in the curvature produce an exponentially 
decreasing correction [IB ] 



V = — (1 - exp(-2or)) , a 2 = 3a + (3 , (5.1) 

r 

backreaction effects due to massless quantum fields induce a power law correction of the 



form |gg7 



y = _GM + a^M^ (52) 



This is the dominant correction at large distances. The numerical constant a depends on 
the matter content of the theory and includes a contribution from gravitons. 

We have a similar picture for the pp-gravitational shock wave. The profile function in 
the quadratic theory is given by 



fan + f Q = f GR {p) - 8G P K (p/^f3) , (5.3) 

and the correction to the GR result is exponentially small at large distances. However, we 
have shown that quantum massless fields induce a correction of the form 

f(p) ~ f GR (p) + ; p^oc, (5.4) 

P 2 

which again is the leading correction at large distances. 

If one treats General Relativity as an effective field theory of an unknown quantum 
theory of gravity, it can be shown that the leading quantum corrections in Eq. ( |5.2| ) will not 
depend on the details of the more fundamental theory: they are a consequence of quantum 
gravity, whatever this theory may be [171]. The same argument applies to the profile function 
of the gravitational shock wave: we expect quantum gravity effects to produce large distance 
corrections proportional to C/p 2 . The coefficient C should be determined by a more careful 
calculation that should take into account all massless fields in nature including graviton 
contributions to the renormalized energy-momentum tensor. However, this contribution 
should be qualitatively similar to the one computed for massless scalar fields in this paper. 

Since we are considering exact solutions to the semiclassical Einstein equations fll.lD , 
it is worth discussing here the question of the runaway solutions that such higher order 



equations may have. In Ref. []I8] it was proposed that a way of getting rid of such unwanted 
runaway solutions is to only consider first order perturbative solutions in the corrections to 
general relativity (here proportional to A). This is clearly too restrictive in the framework 
presented in our paper. In the 1/N approach (N — > oo) the semiclassical Einstein equations 
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hold exactly. In fact, since the scalar fields are taken to be free (i.e. non-self-interacting) 
only the leading order terms in these scalar fields arise. Consequently, one is not only 
allowed to go further than the first order perturbative regime, but one should also take into 
account all the "well-behaved" exact solutions to the semiclassical fields equations. The 
well-behaved solutions we found have been selected from all possible ones by allowing them 
to be represented by a Fourier transform. 
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APPENDIX A: ULTRARELATIVISTIC COSMIC STRING 



It is also of interest to consider the case of an ultrarelativistic particle in 2 + 1 dimensions, 
since its results are also valid for a straight cosmic string boosted in a direction perpendicular 



to it 1141. 

In this case the field equation simply takes the form of Eq. (|3.7Q with replaced by 
dy, where y is the Cartesian coordinate perpendicular to both the direction of motion and 
the axis along the string. 

Consequently, the Fourier decomposition of the profile function leads to 



f(\v\) = f GR (\y\) + ltorGp 



dk 
2^ 



Jk\y\ 





p- 






1 - 


(3k 2 


+ Ak 2 In ( 





(Al) 



where now p represents the momentum per unit length of the cosmic string and Jgr 
—8irp\y\ is the general relativistic result found in Ref. 
By splitting the integral in Eq. (|Al|) we obtain 



f(\y\) = f GR (\y\) + 8Gp\ 



+2VA / dkcos(k\y\*) 
Jo 



r-oo C0s(fc|-y|*) 

1 - \n{k 2 /A* 



;i + fc 2 )(l + fc 2 ln(fc 2 /./l*)) 



(A2) 



where the dimensionless variables have the same meaning as given after Eq. ( |4.4| ). The first 
integral can be found in tables |15|] and its result is ?rexp{— \y\/y/A}/2. Thus giving the 



quadratic theory correction as found in Ref. [13]. The second integral has to be performed 
numerically as we did before for the particle in four dimensions. 
It is easy to check that 



n\y\) 



-IGGp lT 
dk 



A 



cos(fc|y| 



1 + k 2 \n(k 2 /A* 



(A3) 
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is finite at \y\ = 0. The Riemann tensor is then also finite (apart from the 5(u) behav- 
ior). Note that this is already true in the quadratic theory; while in general relativity, the 
divergence was as 

The behavior for large \y\ can be computed on the same lines as before, and we obtain 
f(\y\) ~ -87rGp\y\ + . | y |^oo. (A4) 

\y\ 

Once more, we see that quantum effects due to massless fields are the leading corrections at 
large distances. 
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FIGURES 



FIG. 1. The dashed line represents the path of a null geodesies in the (u,v) plane. Since for 
the spacetime is flat, geodesies are straight lines. At u = they have a discontinuity jump 
described by Eq. ( |2.2| ). The point-like source p is located at the origin and the plane u = contains 
the gravitational shock wave. 



FIG. 2. The 'spatial refraction' of a null geodesies as described by Eq. (|2.3j). Here 9i n and 6. 
are the incident and refracted angles respectively. 



re 



FIG. 3. The higher order corrections to the profile function of the shock wave as defined in 
Eq. (|44]) in units of 8pG. Here p* = p/VA = 27Mp}p and A*=Ap 2 exp(j3/A) are convenient 



dimensionless variables. The correction reaches a finite value at p = and vanishes as 1/p* 2 for 
large p* . 



FIG. 4. d 2 f /dp 2 as defined in Eq. (|4.6| ) over 8pG. This is proportional to the non-vanishing 
Riemann tensor components and show that they diverge as m[ln(p*)] for p* — > 0. 
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